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Abstract. We investigate the asymptotic properties of the large deviation function 
of the integrated particle current in systems, in or out of thermal equilibrium, whose 
dynamics exhibits anomalous diffusion. The physical systems covered by our study 
include mutually repelling particles with a drift, a driven lattice gas displaying a 
continuous nonequilibrium phase transition, and particles diffusing in a anisotropic 
random advective field. It is exemplified how renormalisation group techniques allow 
for a systematic determination of power laws in the corresponding current large 
deviation functions. We show that the latter are governed by known universal scaling 
exponents, specifically, the anomalous dimension of the noise correlators. 



PACS numbers: 05.40.-a, 05.10.Cc, 05.70.Ln, 05.60.-k 

1. Introduction 

1.1. Motivations 

Recently an important concept for the analysis of systems driven out of thermal 
equilibrium, namely that of large deviations in the associated probability distributions, 
has again found considerable attention, both at the experimental and theoretical level. 
Mathematically, large deviation functions have been defined and introduced long ago, 
but only in the early seventies were they shown to play a central role in characterising 
the properties of non-equilibrium steady state (NESS) measures, see [1] for a review. For 
some time this idea remained mostly confined within the mathematical community until 
it was exploited in physics at the theoretical level by Evans, Cohen, and Morriss who 
showed numerically that the temporal large deviation function of a certain observable, 
namely the shear stress, in the NESS of viscous gases under shear displayed a surprisingly 
simple symmetry property [2]. This symmetry property, now known as the fluctuation 
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relation, was then formalised into a theorem by Gallavotti and Cohen [3] . Over the past 
years, a remarkable amount of research has been devoted to the investigation of the 
fluctuation relation in detail, with the goal to explicitly determine the large deviation 
function at the core of the theorem's formulation [U O El [7]. It must be emphasised 
that at about the same period, experimentalists also started to employ, as a tool for 
analysing their data, the concept of large deviations [8j. Their initial interest was 
grounded in the belief that, much in the same way as the intensive free energy per 
degree of freedom constitutes a central quantity if one wishes to relate two systems in 
thermal equilibrium, large deviation functions might similarly allow for a comparison of 
NESS (i?/nam«ca/ properties, yet be applicable to systems both in and out of equilibrium. 
Indeed, it has been a highly desirable, but as yet largely elusive goal of non-equilibrium 
statistical mechanics to be able to classify dynamical systems (at least in a steady state) 
according to a restricted number of (macro- or mesoscopic) criteria. In particular, it 
was expected that large deviation functions of observables such as energy flow, particle 
current (or more generally entropy flow) would display a sufficient degree of universality 
to render this approach fruitful (see refs. |1] through |14j). 

On the theoretical side, a series of recent results however lends support to such 
aspirations. By studying systems of interacting and diffusing particles governed by some 
Markovian dynamics, subjected to an external chemical potential gradient, Bodineau 
and Derrida P,[ini[II], as well as Bertini, De Sole, Gabrielli, Jona-Lasinio, and Landim 
[121 ttSl El; found a general strategy to completely determine the large deviation of 
the particle current flowing through the system. Note, however, that in this setting the 
applied field (and the related particle current) become vanishingly small as the system 
size is increased. Then the only two ingredients entering the final expression for the 
large deviation function are the diffusion constant and the equilibrium compressibility, 
which is an astonishingly simple result: The large deviations of the current are fully 
characterised by properties of the system at equilibrium or in the linear response regime. 

There are two restrictions, however, that apply to this series of recent advances. 
Namely, first they all concern systems exhibiting normal diffusive dynamics. In physical 
terms, this means that the continuum limit can be taken by scaling time according 
to r ~ proportional to a microscopic length scale a squared (often referred to as 
the hydrodynamic limit). There are of course numerous situations in nature where 
diffusion becomes anomalous, which means that the time and length scales appropriate 
for taking the continuum limit are related by r ~ with z ^ 1. The purpose of 
the present work is to investigate the particle current large deviation functions in such 
physical systems governed by anomalous diffusion; specifically, we shall address the 
following models: (i) mutually excluding particles in a driving 'electric' field [l5l[T6]; here 
2; = 3/2 in one space dimension, which is equivalent to the scaling for the noisy Burgers 
equation [T7j, in two dimensions one finds = 2 with logarithmic corrections along 
the drive direction, while normal diffusive behaviour ensues transverse to it, = 2; 
(ii) critical driven diffusive systems (driven model B) fiSl fT9]. for which one similarly 
obtains z\\ = 12/(11 — d) in dimensions d < 5, with again logarithmic corrections 
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at the upper critical dimension dc = 5, and z± = 2; and (iii) particle diffusion in a 
quenched random velocity field, where according to the type of disorder, both super- 
and sub-diffusive behaviour may occur [201 Elj- Based on a field theory representation 
of the non-linear Langevin-type stochastic differential equations representing the above 
processes [22], [231 EH EH] , renormalisation group (RG) methods will prove instrumental 
in extracting the appropriate asymptotic scaling variables in the long-time and large- 
system limit, and will be shown to provide the relevant tool for taking the proper 
continuum hmit [26], [27] . 

The second restriction for most recent investigations is that these studies apply 
only to systems that are weakly out of equilibrium, since, as already mentioned before, 
the applied external field scales as the reciprocal of the system size. To the best of our 
knowledge, it is only for the one-dimensional exclusion process with periodic boundary 
conditions, with [28] or without a bias |7j in particle propagation, that the current 
fluctuations have been accessed analytically to date. The efforts of the present paper 
therefore bear on systems with anomalous diffusion in dimensions d > 1, some of which 
driven far away from thermal equilibrium. 

In the course of this work we shall investigate the properties of the current 
large deviation in three different systems, which we define succinctly in the following 
subsection. In sections 2, 3, and 4 we present the RG analysis and our results for each 
of these model systems. Section 2 itself will also expand on the field-theoretic methods 
upon which we shall rely, whereas most of the technical details will be skipped in the 
subsequent sections 3 and 4. Our conclusions are gathered in section 5. 

1.2. Presentation of the model systems 

We will first investigate the properties of particles (without source nor sink) subjected 
to a driving force along one spatial direction (denoted '||' in the following), whose local 
density p(x, t) evolves according to the continuity equation 9tp(x, t) = —V ■ j(x, t), 
with a particle current j = —D'Vp + pu(p) + rj constructed from phenomenological 
considerations. The first term here is a Fickian diffusion current, the second contribution 
originates from the driving force and describes a density- dependent velocity field u(p) 
along the drive, and finally rj represents Gaussian white noise which accounts for fast 
degrees of freedom not taken into account explicitly in our mesoscopic description. The 
former continuum Langevin description arises for instance from a discrete lattice gas 
description, which we shall use in the following to determine the particular form of the 
velocity field u(p). If we transform from a lattice occupation number representation 
(rij = 0, 1) of the discretised system (with lattice sites labeled by the index i) to 
Ising variables (pi = Ui — 1/2 = =Fl/2, particle-hole {Z2) symmetry dictates that to 
lowest order in the coarse-grained local density field 0(x, t) = p(x, t) — po we have 
M||(0) = e(l — 0^) + . . .. Denoting by jo the average current, renaming e = Dg/2, and 
allowing for a different effective diffusivity XD along the drive direction, we thus arrive 
at the following expressions for the transverse and longitudinal (component along the 
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drive) currents 

j±-j±o= -DV^(P + r)^ , (1) 

Jll-Jlio = -ADV||0-^02 + r^ll . (2) 

The parameter g accounts for the interaction of particles along the driving field. After 
rescaling such that the transverse current sector is effectively in equilibrium, the noise 
correlations read 

{v±^Mr]^,{y,t')) = 2D5,,5^''\^-y)5{t-t') , (3) 

(r/ll(x,t)r^ll(y,t')) = 2Da 5^"^^ - y) 6it - t') . (4) 

Collecting everything into a single stochastic Langevin equation for the density field (f), 
we at last arrive at 

9*0 = D (Vi + A Vj) + ^ V|| 02 + C, (5) 

(C(x, t) C(y, t')) = -2D (Vl + a Vj) 5('^)(x - y) 5(t - t') , (6) 

where ( = —V ■ rj. Further motivations for this phenomenological approach, and a 
detailed study of the scaling properties of this model can be found in [29| [T5| [T6| [T8| [27j . 
We emphasise that in the presence of a drive, and with periodic boundary conditions, 
([5]) and ([6]) in general describe a non-equilibrium system. In one dimension, they reduce 
to the noisy Burgers equation [T7] . 

For our second example, we turn to the standard model of a driven diffusive 
system that exhibits a continuous non-equilibrium phase transition in its stationary state 
[HI [19]. In this case, the transverse current looks like the one for the time-dependent 
Ginzburg-Landau model with conserved order parameter, 

j±-j±o = -I^Vx(r-Vl + ^02)0 + r7^ , (7) 

while the current component along the drive direction is still given by ([2]), and the noise 
correlations by ([3]) and IQ. For the field 0, this yields the Langevin equation 

72 rT2 , ^ J.2^ A 1 \ n tt2 a , ^9 



dt(t> = DVi (^r - Vi + - 02 j + AD Vj + ^ V|| 0^ + C, (8) 

with the noise correlator (Q. The system now has a critical point at r = (within the 
mean-field approximation), which in the absence of the drive g = just describes the 
second-order phase transition in the Ising model with Kawasaki dynamics (model B). 
In the presence of the drive {g 0), however, the universality class changes; in fact, u 
becomes irrelevant for the critical properties of the driven model B, which are entirely 
governed by the non-linear driving term [T8] . 

Finally, we consider non-interacting particles diffusing in a random velocity field, 
as defined by Honkonen [20]. At the mesoscopic level, the system is described over a 
(i-dimensional continuum by a local density fluctuation 0(x, t) subjected to a random 
field '0(xx) directed along ey. One primary feature of the model is that the amplitude of 
the field depends only on the transverse coordinates x^. The density field 0(x, t) evolves 
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according to the Langevin equation dtcf) + V ■ j = with a current j that accounts for 
the random field ^'(xx) solely in the parallel direction: 

j^ = -D^Vx0 + r7^ , (9) 

Jll = V||0-^0 + r7|| . (10) 

The thermal diffusion of particles is encoded by the white noise rj which enters 
the expression for the current. Its correlations ensure that both the transverse and 
longitudinal sectors are in equilibrium, 

(//±,.(x,t)r/^j(x',t')) = 2I^±M(x'-x)5(t'-t) , (11) 
(r/ll(x,t)r7||(x',t')) = 2D||(5(x'-x)5(t'-t) . (12) 
Depending on the spatial correlations of the random field V'(x±) ey, particles can behave 
either sub- or superdiffusively [211 ED] • We restrict our analysis to the case of a spatially 
uncorrelated random field with Gaussian distribution, whose correlations are described 
by a constant amplitude A: 

(V^(x^, t) V;(x'^, t')) = \ 5(x'^ - x^)(5(t' - t) (13) 

This has the consequence that particles behave superdiffusively in the longitudinal 
direction, with a dynamical exponent z\\ = 4/(5 — d) in dimensions d < d^ = ?> [20]. In 
contrast to the previous systems, the superdiffusion here stems from the anomalous, but 
equilibrium wandering of particles in the quenched random velocity field, rather than 
from interactions that drive the system out of equilibrium. 



1.3. Large deviation functions 

In each of these systems the goal will be to determine the probability distribution P{Q, t) 
of the total integrated current in the longitudinal direction Q{t) = /Jdt'/ d'^x j\\{x,t') 
up to time t, and the related large deviation function (LDF) n{q) defined by 

7r{q) = Urn - In P{Q = qt,t) , (14) 

t^oo t 

and expressed in terms of the fiuctuating time-averaged current q. For finite-size systems 
7r(g) is peaked around its maximum qst = limt^oo{Q) /t, which becomes sharper as the 
system size increases (i.e., fiuctuations decrease with system size). In most physical 
systems, 7i{q) does not edgewise reduce to a parabolic form — in other words it describes 
current fiuctuations far beyond a Gaussian distribution, even in equilibrium. As such, 
it appears as a promising tool to probe the distinctive features of NESS as compared to 
equilibrium states. 

For systems driven out of equilibrium by a field E, the fiuctuation relation takes the 
generic form TT{—q) = 7r(g) — Eq. It connects the probabilities of observing longstanding 
current deviations both close to and far from the steady average gst- 

The fiuctuations of the longitudinal current are also fully encoded in the moments 
or the cumulants of Q. We first introduce the generating function 

Z{s,t) = (e-^«W> = ^exp j^t' j d'xj||(x,t') ^ ; (15) 
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moments of Q are then given by the derivatives of t) at s = 0: 
After resorting to the standard Janssen-De Dominicis mapping to a field theory 



= mr) ■ (16) 

s=0 



m [25], we arrive for the systems presented above at the generic form 

Z{s,t) = [ D0P0e-^['^''^'*'*] . (17) 



The specific expressions for the actions S[(j),(f), s;t] will be our starting points in the 
analysis of the systems presented above. We utilise in our work the action functional 
field theory developed in [22l [23], [2l] , rather than the operator formalism of [30] . These 
related but distinct methods lead to quite different techniques, as for instance recently 
exemplified in |31] . 

A few remarks are in order before we proceed with evaluating the large deviation 
function for specific systems. Note that the quantity Z{s,t) defined in f[T5l) . f[T7|) does 
play the role of a partition function, whereas the usual dynamic 'partition function' 
Z{s = 0, t) = 1 (which expresses the conservation of probability) and carries no relevant 
information. Moreover, since the action defined in ([T7|l does not describe a Markov 
process anymore, it is expected that Z{s,t) will grow exponentially with time. It is 
therefore conventional to introduce a time intensive 'dynamical free energy' 

/i(s) = lim - In Z{s,t) , (18) 

t^oo t 

which is related to the current large deviation function 7r(g) as defined in f[T4|) through 
a Legendre transform, 

Tc{q) = maxs{fi{s) + sq} . (19) 

Formally, the dynamical free energy fi{s) is nothing but the generating function for the 
cumulants of the integrated longitudinal current Q in the asymptotic long-time limit: 

=.taM):>.. (20) 

ds" s=0 t^oo t 
1.4- Steady states with non-zero value of s 

In analogy with standard thermodynamics, we observe that the quantity Z{s,t) 
introduced in fllSp facilitates the description of current fluctuations in a canonical way 
(i.e., at fixed intensive parameter s) rather than microcanonically (at fixed value of the 
time-extensive quantity Q). 

Besides constituting a useful computational trick, introducing Z{s,t) also provides 
us with a physical picture of the very configurations that give birth to large deviations. 
Although at s 7^ the action defined in f[T7|l does not correspond to a stochastic 
process anymore, one can understand S[(l), 0, s; t] as describing not a single, but rather 
an ensemble of systems evolving in parallel, with dynamical rules that favour histories 
with non-zero longstanding deviation of the particle current [32]. The value of physical 
observables such as correlation functions inferred from S[(j), (j), s; t] at s 7^ can thus be 
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understood as the typical value of these observables in a modified steady state wherein 
the mean current is enforced to take an average value q{s) that is different from the 
average q^t taken in the steady state. We see for instance from ( |T5l) that negative values 
of s favour histories with an excess current in the direction of the field. Quantitatively, 
the correspondence between q[s) and s is the same as in the Legendre transform ( JT9l) . 
The reader is referred to [32| [33] for further examples and for a numerical exploitation 
of this picture. 



2. Driven diffusive system with mutually exclusive particles and the noisy 
Burgers equation 

2.1. Field-theoretic formulation 

We wish to determine the distribution function of the integrated longitudinal current. 
As outlined in section [L3l this amounts to determine the 'dynamical free-energy' Z(s, t) 
defined in (JTSll . rewritten as Z(s,t) = J T)(f)T) (j) Q-^i'l'''l'''^'A^ Using periodic boundary 
conditions we find that the dynamical action S[(j), (p, s; t] is given by 

S[(j),(j),s;t]= j (i'^x j dt' (f){dt' - DVl- DXVl)(j)- D {V±(j)f - Da 

-^0' + ^(V||0)0^] +sjoLH-s'DaL't. (21) 

(We use the Ito convention as regards time discretisation.) It mainly differs from the 
stochastic action S[4>, 0, s = 0; t] through a new quadratic term in the field 0, aside from 
the two deterministic contributions ~ L'^t. 

In the case of non-interacting particles, each of the cumulants {Q^)c is extensive in 
the system volume, which means that the corresponding free energy /iindcp(s) takes the 
scaling form 

/i,„,.,(s, L) = fi,^,^^{s) = {B^s + B^s'' + . . .) , (22) 

where /iindep('S) has a well-defined power expansion around s = 0. It is known from 
specific examples, in |35l [7] or out of equilibrium [28],|M|, that interactions can lead to 
an expansion of /i(s) = L~'^fi{s,L) involving non-integer powers of |s| (in the infinite 
volume limit). Such non-analytic behaviour at s = simply mirrors the existence of 
an infinite cumulant of the current as L — > oo. For instance, the LDF for the totally 
asymmetric exclusion process (TASEP) in the large system limit (with finite density p) 
reads [28l El] 

/iTASEp(s) = -p(l-p)s + ^4^[p(l-p)]t|s|t + ... fors<0, (23) 

5 

which means that the variance of the particle current grows faster than the system size, 
whereas in the symmetric exclusion process (SEP) one has [351 C] 

/isEp(s) = ^p(l - p)s' + ^P^^ -p)]^s\I + ... for s < . (24) 
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The anomalous scaling of the current LDF occurs here for a higher-order cumulant, 
reflecting that the anomalous distribution of the current arises from sharper details 
than in the TASEP. From a general point a view, we expect the exponents in ( l23l) or 
(!24|) to be universal (for other models within the same universality class), while the 
prefactors should be model-dependent. 

In our case, the continuum and large system size limits will be tackled in section 
12.21 with the presentation of the - unavoidably technical - RG analysis of the driven 
diffusive system we are considering. We then proceed to compute the dynamical free 
energy in subsection 12.31 before extracting its asymptotic behavior as a function of 
s inEl 

2.2. Renormalisation 

The presence of the additional quadratic vertex in S'[0, 0, s; t], in spite of allowing for new 
one-loop diagrams contributing to the three-point vertex function r*^^'^\ does not alter 
the remarkable property of this model that there are strictly no singular loop corrections 
to r^^'^), i.e., all new contributions lead to singularities that cancel out precisely. One 
can explicitly check that order by order in the perturbation expansion, or infer this 
property from Galilean invariance, see below. The same statements hold for r(°'^\ 
which receives no singular loop correction either. However, both F*-^'^) and F*^^''') have 
vertex corrections given by exactly the same diagrams as with s = 0, but with values 
that are now s-dependent. As a consequence, no new renormalisation is required as 
compared with the stochastic case s = 0. Note that for cr = A the Einstein relation also 
holds in the longitudinal sector. Then the symmetry 

= S[^{-t) - , (25) 

which for s = arises from the existence of a free energy functional wherefrom the 
Langevin equation derived, continues to hold for s 7^ - in spite of the action not 
representing a stochastic process anymore - with the consequence that the vertex 
functions F^^'^^ and F^^'^^ are given by their tree-level expressions. 

In the remainder of this subsection, we recall for later use the main lines of the 
renormalisation analysis, which is the same as in the original approach [151 [27]. Denoting 
by K an arbitrary momentum scale, the scaling dimensions of the various fields and 
coupling constants are fixed by the requirement and convenient choice that the action 
S as well as the diffusivities D and D\ be dimensionless, 

[x] = K-\ [t]=K-\ [q] = K, M=k\ [<P\ = [$\=k'I\ (26) 
[D] = [a] = [\] = n\ [g]=K'~'^l\ [s] = k^+'I^ (27) 

Of course, any other choice for the scaling dimensions of the field and the diffusion 
constant D would eventually lead to the same physical consequences. We infer from 
the scaling dimension of the nonlinearity the upper critical dimension = 2, and also 
[s g D] = K^. Next we define renormalised fields and parameters via 

0H = <V, ^n = Zf^, (28) 
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propagator: 



noise vortex: 



nonlinear vertex: 




Figure 1. Propagator, noise vertex, and 
expansion for the driven lattice gas. 



1 

-iLt;+D(q^+Ayp 

2D(qi + agp 



\gD igy 



nonlinear vertex used in the perturbation 



D^ = Zr>D, a^ = Z,a, X^ = ZxX, g^ = Zgg k-^+^'K (29) 

In the transverse sector, i.e. on the q\\ = momentum shell, the action is merely 
Gaussian and we have to all orders in the perturbation expansion 

r^'''\q^,qil=0,uj) = tuj + Dql, (30) 
r(''°Hq±,?||=0,a;) = -2Dqi. (31) 

As a consequence, Z^^'^ Z^^'^ = 1 and Z^ = 1. Then, exploiting that the transverse 
sector is in equilibrium, we obtain from the associated fluctuation-dissipation theorem 
the relation Z^ = Z^^^'^ Z^^'^. Thus we see that neither the fields nor D renormalise: 

= = Zd = 1 , (32) 

and consequently infer the exact transverse Gaussian scaling exponents 

r] = , z^ = 2 . (33) 

The invariance of the action under Galilean transformations 

0'(x'_i_, t') = 0(xx, — Dgut, t) — u (34) 

shows that cj) and u must transform in the same way under renormalisation {u is the 
constant velocity of the Galilean transformation). Moreover, the product Dgu also 
remains invariant under RG. Thus, 

Z, = Z^' Z^'/' = 1 , (35) 

and the coupling constant g is not renormalised either. 

We now have to perform an explicit loop expansion using the propagator and 
vertices depicted in Fig. [H We use the dimensional regularisation scheme and perform 
the expansion in powers of e = dc — d = 2 — d. A straightforward explicit computation 
of the two-point vertex functions yields 

Zx = l + ^i3 + w^) + Oivl) , (36) 
Z^ = l + ^ {3w~' + 2 + 3w^) + Oivi) , (37) 
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where we have defined renormahsed effective couphngs 

73/2 g"^ ^ d-2 n r(2-rf/2) 

Wh = 3^, v^ = Z^' ^CdK , with Crf = ^^3^-;^ , (38) 

and Wilson's flow functions 

7a = /^^kIo In ^ = (3 + iUr) , (39) 
A lb 

7. = «:9,|o In — = (3<i + 2 + 3w^) (40) 

to flrst order in t>R (the subscript '0' here indicates that the derivatives are to be 
performed in the unrenormahsed theory). The RG beta functions follow 

Pw = fid^lo u7r, = Wj, (7^ - 7a) = {w^ - 1) (u^R - 3) , (41) 

Pv = K(9k|o t^R = v^(^d-2-^'j^ . (42) 

RG flxed points are now determined by the zeros of these beta functions. At any non- 
trivial flxed point < i;* < 00, ty* = 1 is stable and the system reaches effective 
equilibrium (cr* = A*). To all orders in perturbation theory we then obtain from ( H2l) 

ll = i: = lid-2) , (43) 

in dimensions d < dc = 2. 

The dynamical correlation function scales as 

C{c,,.) = ciT--'^^c(j-^„^] , (44) 



|q±r"" |q±i 



A = = . (45) 



where in < 2 dimensions exactly 
7i^ 2- 
2 3 



We thus obtain 

and as expected, zq) = 3/2 in one dimension. For d > 2 the system scales towards the 
Gaussian flxed point Vq = 0, whence one arrives at the mean-fleld scaling exponents 
A = 0, 211 = = 2. 



2.3. Evaluation of the cumulant generating function fi{s) 

Having recalled the main results of the RG analysis, we can now turn to the evaluation 
of the dynamical free energy /i(s). The dynamical action (pTj) contains an extensive (in 
time as well) deterministic contribution, L'^t [sjo — s'^Da), and a fluctuating one, for 
which it will prove sufficient to extract the tree level approximation. Quite unfamiliarly 
in dynamical field theory, the latter contribution arises from the normalisation of the 
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path integral. The Gaussian contribution to the action S'[0, 0, s; t] takes the form 



with, conveniently written in Fourier space, 



ro(q,tu) 



-2D {(\\ + aqfj -iuo + d{(\\ + Ag[ 



icu + D qj_ + \q\ 



'Dgs 



We have 



^^^(q) =D^c^l + XqiY-2sgD^cil + aql) , 
and integrating over the fields yields the following tree level contribution to fi{s) 

1 



D(j)'D(j) exp — 



exp 



exp 



exp 




ro(q,u;) 



ro(q,u;)|,=o 
2sgD\cil 



2 + D2(q2^ + Ag2)2 



Ag;)'-2.^?(qi + agj)-(ql + AgJ) j 



(47) 

(48) 
(49) 

(50) 
(51) 
(52) 

(53) 



At the RG fixed point the values of A and a are not fixed, but become equal 
[w* = 1). Without loss of generality, one can thus start from equal bare coupling 
constants cr = A (as was for instance implicitly assumed in [29]). The last integral ( l53ll 
can now be carried out explicitly, and its infrared divergence is isolated by extracting 
the first term of the small s expansion. This gives 

ri . . 2 1 D\-^/^(-gsy+'^/^T(d)T(-l-d/2) , , 



2d(27r)'^/2 [r(rf/2)]^ 



2.J^. Scaling behaviour of fi{s) 



To determine the expansion of /i(s) in powers of s in the continuum limit, and at the 
stable RG fixed point = 1 (A = cr) we notice that making explicit the 1/e pole in 

C,g' i-gs)-^/'- 



1 + 



(55) 



1 + ^ and finds (at the convenient normalisation 



A3/2 4£ 

one recognises the expression for Zx 
point K = (—sg)^) 

L-''fi{s) = -joS + DX^s\ (56) 

Finally, using that in the vicinity of the RG fixed point A^ ~ Ak'^J , with the anomalous 
scaling dimension 7^ = 7* = — |e of the noise correlator, the asymptotic scaling of /i(s) 
as a function of s reads (below the critical dimension dc = 2) 

L-V(s) = -Jos + A,s'-^/' = -Jos + A^s^'^'^/' . 



(57) 
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The constant prefactor Ad is non-universal, contrary to the exponent 2— e/S = ((i+4)/3, 
which holds exactly to all orders in perturbation theory. In one dimension, we recover 
the exponent 5/3 of the totally asymmetric exclusion process [2HI IM], see ( 123|) . This 
results expectedly corroborates that the TASEP and systems described by the Burgers 
equation belong to the same (dynamical) universality class. 

We emphasize that the scaling exponent in (|57|) merely follows from the noise 
renormalisation and could thus have been obtained through straightforward dynamic 
scaling. A simpler, though less explicit, way to arrive at that conclusion is to employ 
the matching condition stemming from the A'' IR divergence in (15^ : The only factor 
acquiring an anomalous dimension in the crucial second term in (15^ is the parameter 
cr; with 7* = | (c? — 2) and the matching condition k ~ {—sgY/"^ that directly follows 
from the scaling dimensions, we again arrive at fi(s) + jos ~ s2+(<i-2)/3_ Jn^jgg^j^ ^^j^jg 
will be our method of choice in the following subsection as well as for the remaining two 
systems. 



2.5. Logarithmic corrections to fi{s) in two dimensions 

Precisely at the upper critical dimension dc = 2, the RG flow equations provide access 
to the logarithmic corrections to the mean- field critical power laws. In our context, they 
also specify the logarithmic correction to the term in fi{s). To determine how the 
renormalised coupling constants vary under scale transformation k ^— > in, let us recall 
that running couplings A(£) and v{i) are defined through 

«^ = o.(^)Am. ^^=Am. (58) 

with initial conditions A(l) = Ar, v{1) = and where 7a(^) = 1\{v{i)), [3^{t) = 
Pv{v{i)). In dc = 2 dimensions, the flow equation for v{£) reads 



which is solved by 



-ij'{i) + Oiv'ii)) , (59) 

(60) 



1 - In 



Upon inserting this result in the flow equation for X{i), one obtains 

A(£) ~ A(ln£)'/' , (61) 

which gives the precise form of (!57|) at dc = 2: 

L'' /i(s) = -Jos + A2 s\ - In \s\f' . (62) 

Of course, such logarithmic corrections are familiar within the framework of RG 
calculations at the upper critical dimension. Yet we stress that this result constitutes the 
first exact result for the LDF scaling of a two-dimensional system, where the methods 
of integrable systems fail to apply. 
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Figure 2. One-loop graphs which renormalise the correlation function in the s-state. 



2. 6. Correlation functions 

In this subsection, we aim at determining how the correlation function are modified 
for small (negative) values of s. As announced in section 11.41 states with s < 
are characteristic of histories with excess current in the direction of the driving field. 
Let G^™'") = (0^0") denote the m-0 and n-0 correlation function. For the complete 
quadratic part of the action, we find that 

causal part non causal part 
. ^ s . ^ s 



G(^'i) (q, t, t') = e(t' - 1) ^ e-^^(*'-*) + e(t - t') 



s 



with the abbreviation fH^ . and with the convention that the step function satisfies 
0(0) = 0. (Note that the above calculations hold only for s < 0). The diagrammatic 
expansion contains no causality constraint on the loops anymore. The long-range spatial 
correlations in the s-state are obtained from an expansion in q at minimal order. To 
one loop, the correlation function is renormalised through the graphs shown in Fig. [2j 
To lowest order in q, we find that all contributions cancel except for the single loop 
diagram depicted in Fig. [3l which yields 

C(q, t, t) - C(q, t, t)l^, = {-2gs)-'l' (qi + Agf ) 

l + /A-i/^(-^s)-/^ 



X 



.-el2 1 ^11 



(64) 



IQixe q5_ + Agji 

Averaging along the equilibrium directions, i.e., at = 0, one recognises 

C(g|| , 0, t, t) - C(g|| , 0, t, t) = {-2gs)-^'^ X'J' g,, , (65) 
and finally obtains the exponent of s in the s-correlation function: 

C(g|| , 0, t, t) - C(g|| , 0, t, t) (X i-s)-'/'-^/' . (66) 
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cm 




Figure 3. Single diagram contributing to C(q, i,t) in the s-state, at minimal order 
with respect to the external momentum q. The expressions for the propagators are 
given in 



In one dimension, this gives the power law (— s)^^/^, which, to our knowledge, has not 
been obtained through exact approaches for systems in the Burgers/KPZ universality 
class. 

3. A driven diffusive system with a continuous phase transition 

3.1. Renormalisation 

The previous considerations and calculations can readily be generalised to the critical 
driven diffusive system described by the nonlinear Langevin equation ([8]) with noise 
correlator (jS]). Since only the transverse sector in momentum space is rendered 
critical, we merely need to replace the propagator in Fig. [1] with the expression 
l/[—iuj + £)q^(r + q^) + D\ gjj]. Dimensional analysis yields [r] = /t^, and at criticality 
(r = 0) [g||] = [q'x]^ = K^, and hence [uj] = k*^. This yields for the scaling dimension of 
the nonlinear coupling [g] = k^^"'')/^, wherefrom we infer the upper critical dimension 
dc = 5. 

Renormalising this model proceeds in much the same way as described above in 
section [22] (see [iHl ET]). Since again the nonlinear vertex is proportional to iq\\, the 
transverse sector is Gaussian with 

r^'-'^ (qx, QW =0,uj) = iuj + Dql(r + ql) , (67) 
T^'''Hci^,qil=0,uj) = -2Dql. (68) 

Consequently ( !32l) holds here as well, implying that 

V = 0, iy=^, z^ = A (69) 

exactly. The system is of course still invariant with respect to Galilean transformations, 
whence (1351) is valid. The RG beta function for the effective nonlinearity v = g^j}?!'^ 
now reads 

ii. = vAd-h-\-iy\ , (70) 



2 

and thus again to all orders in the perturbation expansion 

ll = lid-5), (71) 
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and consequently 

3.2. Evaluation and scaling behaviour of fj,{s) 

The Gaussian contribution to the action S[4>, (p, s; t] now becomes 



exp 



2 J (27r)<^ 27r V + D^[(i\{r + c^^) + \q^^^ 



(73) 



Thus we observe that we require the scahng dimension of the noise strength cr, which 
is an irrelevant parameter, in the RG sense, for this model. A straightforward one-loop 
calculation gives 

7A = -^H, 7^ = 2-^, (74) 

so for e = 5 — d > 0, 

T: = 2 - I + 0{e') . (75) 

However time now scales according to \t\ — k"^, and the appropriate matching condition 
becomes k ~ {—sgY^'^. Thus we finally arrive at 

L-^,x{s)^-jos + Bas'-^l'\ (76) 

with a non-universal amplitude Bd- Note that this result holds only to first order in the 
dimensional expansion for (i = 5 — e < 5. 

At the critical dimension = 5, 

wherefrom we obtain the logarithmic correction 

M = -Jos + i35 s\ - In Islf" . (78) 

4. Superdiffusion in a random velocity field 
4.1. Mesoscopic formulation 

Before embarking on the determination of //(s) for the case of (super-)diffusive particles 
subject to a random velocity field, we have to establish how the random field ip intervenes 
in the definition of the LDF of interest. The physical features of the current fiuctuations 

are contained in the disorder-averaged cumulants {Q{t)"')c- Here and below the averaging 
over thermal disorder will be noted by (...), while the average over the random field will 
be indicated with an overbar, .... We are thus ultimately interested in determining 
the LDF 



//(s) = lim - In Z^{s,t) . (79) 

t—*oo t 



Current distribution in systems with anomalous diffusion: RG approach 



16 



The disorder-dependent partition function Z^{s,t) reads 



Z4s,t) = {exp[-s dt'd'^xj||(x,0]) , (80) 





and as in the previous models, a field-theoretic reformulation of (l9ll-( |TT]) and (IHOl) enables 
us to rewrite the partition function as 

Z4s,t) = J P^P(/.e-^*[^''^'''*] , (81) 

where the quenched and s-dependent action is given by 

5^[0,0,s;t]= -s^D\^LH + J d'^x dt' ^{df - D^Vl- D\\Vl-tljV\i)(j) 

- D^iV^W - D\\{V\\W - #0] ■ (82) 

The form of this action is related to Honkonen's original (s = 0) action [20] through 
the canonical transformation 

= P, = P + P, (83) 
whereupon takes the form 

S^[p,p,s;t] = -s^D\\LH 

+ j ^'^^ dt' p{dti - D^V\- D{^\-ijV{) p + s{p + p)i) . (84) 

For s = this is precisely the action that was studied by Honkonen |20j . 

In a microscopic lattice gas formulation (independent particles diffusing in a random 
field), Z^{s,t) can also be obtained from the standard Doi-Peliti approach [35]. Up to 
irrelevant terms, one arrives again at the action fl82l) . with the field being related to the 
original operators through a Cole-Hopf transformation. This illustrates that, either in 
the Janssen-De Dominicis or in the Doi-Peliti manner, the well-tried functional methods 
of statistical physics can be readily extended to tackle large deviation functions and the 
associated s-modified states, as previously proposed in [57] . 

4-2. Renormalisation 

The theory has two superficially divergent graphs which lead to the same renormalisation 
of D\\. As can be seen from direct inspection of possible diagrams, the additional sip(j) 
term in the action does not lead to further renormalisation, as compared to the s = 
case. We define as in [20] renormalised couplings through D\\ = ZD\\^^, A = ARK^ The 
perturbative expansion is organised in powers of a single dimensionless parameter 
that we fix to 

Examination of the one-loop contribution to F'-^'^-' leads to 

Z = 1 - ^ , I3u = u^{u^ - e) . (86) 
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Below the critical dimension rf^ = 3, the RG flow has a single non-trivial fixed point 
= e, at which the system exhibits superdiffusion [20], namely 

(xj) ~ ; (87) 

in other words, the dynamical exponent in the longitudinal direction is 

4 



^11 



2 + e 



4-3. Determination of ^{s) 

From a dynamical point of view, the study of disordered systems does not pose any 
particular problem since conservation of probability ensures that the dynamical partition 
function remains equal to 1 as time evolves. This has the well-known consequence that 
quenched disorder can be treated without resorting to the replica trick. However, in the 
present case, our interest aims precisely at the nonzero and s-dependent dynamical 
partition function Z^{s,t), which grows exponentially with time, with a disorder- 
dependent rate. The dynamical free energy we pursue, denoted by /i(s), is obtained 
as the average of this rate, see fl79l) . The replica trick represents one convenient means 
to achieve this disorder averaging and to determine /i(s) via 

1 d 



Li(s) = lim 



[Z^{s,t)Y . (89) 

n=0 



In contrast with previous work [38j using replicas to probe dynamical aspects 
of disordered systems, our approach leads to non-trivial couplings between replicas. 
Integration over if) yields 



[Zi.{s, t)Y = J V(j)a e-^"[^'*'^^*] , (90) 
with an effective action Sn[(pa, (pa, s;t] whose quadratic part reads 

^n[0a,0a,s;t]= - nL'^tDys^ - J 0,(q^, 0, 0) 0fc(-q^, 0, 0) 



d'^q duj 



+ y [M-^^ + D ■ q^)0. - (D . q^)0fj , (91) 

where we have set D ■ q^ = -D±(q^)^ + -Dy (9|j)^ for clarity. Integration over the noise has 
given birth to a uniform coupling between replicas of the field 0, in the u = 0, q± = 
sector. However, no new renormalisation is needed, as can be directly checked. 

The free energy /i(s) is obtained as the ratio of the determinant of the action at 
s 7^ and s = 0. In our case, the effective action (I9T1) is diagonal in Fourier space, and 
the only modes which depend on s are those with gy = and u = 0. We thus require 
the matrix determinant 

A„(q J - det I ^ I , ( J2j 
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where I„ denotes the n x n identity matrix, and JI„ represents the matrix whose entries 
are all 1. For later convenience, we have also added a mass m to the (fxp term. Using 
standard results about block determinants, we find 



An(qx) = (-1)" (^±qi) [Dl (qi + mY- 2nXs'D^ql\ . (93) 
The corresponding (tree-level) contribution to fi{s) is then given by 



d'-'q^ An(qx) ^ f d'-\^ ql 



(27r)^-i A„(q^)|,=o D^J (27r)^-i (g2 + ^2) 



,^As^(2-e)m "„/e\„/^ e 



2 an,_ _ 

Combining this result with the deterministic contribution to /i(s) we recognise 
4-4- Scaling behaviour of fi{s) 

As in the previous examples, we infer from fl95l) that the continuum limit in /i(s) is 
determined by the renormalisation of the noise coupling constant. Using j"^^^ = — e, one 
gets 

L-V(s) ~ s^K"' , (96) 

but from the explicit expression 0921) . the normalisation point depends on s only through 
As^. Since [As^] = k^, we obtain the scaling behaviour, valid below the critical dimension 

dc = 3, 

L-V(s) ~ . (97) 

We note that although this system is in equilibrium, the expansion (!97|) implies that 
that in a finite-size system the variance of the current grows with the linear size L as 
for < 3, in contrast to the L'^ behaviour observed in normally diffusive systems. 
Again, this anomalous behaviour is exemplified here for the first time. 

At the critical dimension dc = 3, the running coupling 

D\\{i) ^ D^lni , (98) 

which gives the logarithmic correction 

/i(s) ~ s2(-ln|s|) . (99) 

At the upper critical dimension too, the RG analysis leads to an asymptotically exact 
result. 
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4-. 5. Correlation function 

As well as allowing access to free energies, the replica trick can also be used to determine 
the correlation functions, averaged over thermal fluctuations and the random field. This 
is achieved by noting that 

(0(q^,O,O)0(-q^,O,O))e= (0a(qx,O,O)0,(-q^,O,O)),|„^o . (100) 

Thus we obtain at the tree level the correlation function in the s-state as 

(</)(q^, 0, 0)</)(-q^, 0, 0))e = (l - ^s^-^^ ' (^O^) 

The spatio-temporal correlation function reads 

CM = C(x,i)U— y p^^e-- , (102) 

i.e., the supplementary correlations in the s-steady state are still uniform in time and 
along the longitudinal direction, which means that states conveying a (slightly) atypical 
value of the current do not break theses symmetries of the steady state. Including loop 
corrections to the correlation function does not alter this result, but affects the power 
in (11021) . as a consequence of superdiffusivity. 

5. Conclusions 

In the previous sections, we have shown that the implementation of dynamical RG 
techniques to the already well-studied (see refs [2] through [H]) dynamical free energy 
fi{s) associated with the particle current provides new insights into the singular 
dynamical macroscopic limit of large systems. As is the case for the static free energies 
in thermal equilibrium setups, fi{s) exhibits singularities only in the infinite volume 
limit. In the systems presented in this work, /i(s) becomes non-analytic at s = 0, where 
s is the parameter canonically conjugated to the particle current: 

L-V(s)+joS ~ , with l<a<2. (103) 

In finite-size systems, the variance of the integrated current grows with the system size 
faster than linearly. We have shown for the systems under consideration here that the 
exponents a were simply related to the respective dynamical exponents z through the 
anomalous exponent for the noise strengths. In our examples, the divergence of the 
variance (a < 2) emerges as a direct consequence of superdiffusion {z < 2). 

Our second system provides an example for a current LDF in presence of twofold, 
namely static and dynamic criticality. Our final model system constitutes the first 
example of an equilibrium system for which the fluctuations of the particle current were 
shown to produce anomalous scaling with system size. 

In previous works, which remained confined to a specific one-dimensional system, 
namely the asymmetric exclusion process, an exponent related to a was determined 
using the Bethe ansatz [28l [Ml [33]. Of course, some of these existing results can 
be recovered through our field-theoretic methods employed here (such as the scaling 
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exponent m. d = 1), but we have obtained additional asymptotically exact results in 
higher dimensions. As well as opening the path to a wide range of applications, our 
method also connects to other multi-purpose instruments of statistical physics (such 
as replicas and supersymmetry for disordered systems) that can be deftly adapted to 
describe the s-dependent states. As illustrated in sections 12.61 and 14.51 our method 
also gives access to the correlation functions in the s-state, which to our knowledge 
cannot be obtained from Bethe ansatz computations, but still provides valuable physical 
information about the critical behaviour in the vicinity of s = 0. 

Finally we remark that for the driven systems investigated here we have not been 
able to access the s > states which probe realisations of the process with a current 
opposite to the average current. This is likely because instead of remaining homogeneous 
in time and space, the typical trajectories develop heterogeneities in this situation. We 
may speculate that, as discussed by Bodineau and Derrida [ini [H], non-stationary 
profiles, perhaps with shock-like spatial structures, will play an important role. This 
issue would certainly be worthwhile further detailed investigations. 
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